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ABSTRACT 

The recent model of quasi-periodic oscillations in neutron star binaries 
(Osherovich and Titarchuk 1999, Titarchuk and Osherovich 1999) has suggested 
the existence of two branches of QPOs due to the influence of Coriolis force on 
the linear Keplerian oscillator: one branch with frequencies ~ 400 — 1200 Hz and 
another branch with frequencies an order of magnitude lower. The frequencies of 
the high branch Uh hold a hybrid frequency relation with the Keplerian frequency 
'^K- Vh' — + f^^/TT^, where Q is the rotational frequency of the star's 
magnetosphere. The frequency of the low branch is vl = iS^/T^){.^Kl^h) sin 5, 
where 5 is a small angle between SI and the vector normal to the plane of 
Keplerian oscillations. The observations of the source 4U 1702-42 (Markwardt et 
al 1999) have shown that the centroid of the ~ 35 Hz QPO tracks the frequency 
of the kilohertz oscillations. We interpret the ~ 35 Hz oscillations as ul and 
find 5 = 3.9° ± 0.2°. Our results make 4U 1702-42 the second source (after Sco 
X-1) for which the theoretically derived lower branch is identified (within our 
model) and 5 is calculated. The inferred angle 5 stays approximately the same 
over the significant range of vk (650 - 900 Hz), as expected from the model. 
Based on our model we present a classification of QPO frequencies in the source 
4U 1702-42 observed above and below vl. 

Subject headings: X-rays: bursts-stars: accretion disks — radiation mechanisms: 
neutron — stars:individual: 4U 1702-42, Scorpius X-1 
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1. Introduction 

The discovery of kilohertz quasiperiodic oscillations (QPOs) in the low mass X-ray 
neutron star binaries (Strohmayer et al. 1996; Van der Klis et al. 1996) was followed by 
similar results for nineteen sources (Van der Klis et al. 1998, Strohmayer, Swank and 
Zhang 1998). For most of them, the Rossi X-Ray Timing Explorer (RXTE) mission showed 
the persistence of twin peaks in the spectrum (~ 400 — 1200 Hz). In the lower part of the 
spectrum of Sco X-1, van der Klis et al. (1997) found two branches ~ 45 and 90 Hz which 
correlate with the frequencies of the above-mentioned twin peaks. The nature of the twin 
peaks has been discussed intensively. The apparent constancy of the difference Ai/ for the 
twin peaks in some sources led to the beat-frequency interpretation (a concept originated 
by Alpar and Shaham 1985). Below, we describe difficulties of this model following the 
recent assessment by Mendez and van der Klis (1999). Within the beat-frequency model, 
the higher peak in the kHz range is identified with the Keplerian frequency (van der Klis et 
al. 1996). The lower twin peak in this model occurs as a result of beating between the high 
peak and the spin frequency of the neutron star Vspin which is believed to be observed during 
type 1 X-ray bursts as v^yurst (or half that) (Strohmayer et al. 1996, Miller, Lamb & Psaltis 
1998). Within the original beat-frequency model, Ai/ and Vspin should remain constant. 
The observed 20 - 30% variation of Ai/ (Van der Khs et al. 1997, Mendez et al. 1998, Ford 
et al. 1998, Mendez & van der Khs 1999) undermined the beat-frequency interpretation. 
Convincingly, Mendez and van der Klis (1999) showed that for atoll source 4U 1728-34, Ai/ 
is not equal to Vhurst even for the lowest inferred mass accretion rate. In their words, this 
seems to rule out the simple beat-frequency interpretation of the kHz QPOs in LMXBs and 
some of the modifications introduced to explain the results of a variable Ai/. 

A completely different paradigm has been proposed by Osherovich and Titarchuk 
(1999) and by Titarchuk and Osherovich (1999) (hereafter OT99 and T099, respectively) 
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to explain QPOs in neutron stars. The new model is based on the idea that the twin peaks 
occur as a result of the Coriolis force imposed on a one-dimensional Keplerian oscillator. 
The second section of our Letter contains the assumptions and predictions of the model. 
Verification of the predictions concerning the low frequency branch for the source 4U 
1702-42 is presented in section 3. Discussion and summary follow in the last section. 

2. Oscillations in the Disk and the Magnetosphere Surrounding the Neutron 

Star 

The following assumptions constitute our model (OT99 and T099): 

a) The lower frequency of the twin kHz QPO is the Keplerian frequency 



where G is the gravitational constant, M is the mass of the neutron star and R is the radius 
of the corresponding Keplerian orbit. 

b) There is a boundary layer: transition region between the NS surface and the first 
Keplerian orbit. Within this layer, the radial viscous oscillations are maintained. There 
is also a diffusive propagation process of the perturbations generated in this layer and 
characterized by a break frequency i/b related to Uy. Outside of this layer, the disk has 
Keplerian radial oscillations with (see also Titarchuk, Lapidus & Mushmov 1998). 

c) The magnetosphere which surrounds the neutron star has a rotational frequency ft 
which is not perpendicular to the disk. 

d) inhomogeneities (hot blobs) thrown into the vicinity of the disk participate in 
the radial Keplerian oscillations but simultaneously are subjected to the Coriolis force 
associated with the differential rotation of the magnetosphere. 




(1) 
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The consequences of the model are straightforward. The problem of a linear oscillator 
in the frame of reference rotating with rotational frequency fl is known to have an exact 
solution describing two branches of oscillations 

''^ = 4+(")' (2) 

ul = i^/Ti) (i^K/i^h) sin 6 (3) 

where 6 is the angle between CI and the plane of the Keplerian oscillations. Formulas (2) 
and (3) are derived under the assumption of small 6 (OT99 and references therein). Thus, 
the high frequency of the observed twin kHz peaks has been interpreted as the upper hybrid 
frequency i/h or the frequency of the high branch of the Keplerian oscillator under the 
influence of the Coriolis force associated with magnetospheric rotation. The existence of 
the lower branch (i/^) is a prediction of the model. Oscillations ~ 45 Hz and 90 Hz in Sco 
X-1 reported by van der Khs et al. (1997) are found to flt formula (3) for 5 = 5.5 ± 0.5° 
when interpreted as 1st and 2nd harmonics of vl (OT99, T099). In the low part of QPO 
spectra according to T099 there are two more frequencies (namely the break frequency 
and frequency of viscous oscillations v^) related to each other 

i/fc = 0.041i/^-^\ (4) 

We now consider the application of our model to source 4U 1702-42 studied in detail by 
Markwardt, Strohmayer & Swank (1999). 

3. Determination of 5 for the Source 4U 1702-42 

For Sco X-1, we have shown that the solid body rotation Q, — Q,q — const is a 
reasonable flrst order approximation (OT99). Theoretically, 0. depends on the magnetic 
structure of the neutron star's magnetosphere. Within the second order approximation, the 
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Q profile has a slow variation as a function of uk, which we modeled for Sco X-1 and for the 
source 4U 1608-52 in OT99. For the source 4U 1702-42, there are not enough simultaneous 
measurements of Vh and uk to reconstruct the fl profile. Thus, we restrict ourselves to the 
first order approximation and adopt fi/27r = Qq/Stt = 380 ± 7 Hz, found in OT99. The 
observed and used in OT99 for the calculation of Qo/2n according to formula (2) are 
presented in Table 1 which contains the data kindly provided by C. Markwardt. Then from 
formulas (2) and (3) we have the expression for S, namely 



where uk is the observed Keplerian frequency (the second from the top kHz QPO) and i/l 
is the alleged frequency of the low Keplerian branch. Within the first order approximation, 
the angle S should stay the same for all uk, since 6 is effectively the angle between the 
equatorial plane of the disk and the plane of the magnetic equator. The first two columns 
of Table 2 present frequencies observed by Markwardt et al. (1999) for the source 4U 
1702-42 with the corresponding date of the measurements. These two columns repeat the 
first two columns of the table in Markwardt et al. (1999) with the addition of data for 
July 30, 1997 taken from the text of the same paper. In the third column of our Table 2, 
we give our theoretical interpretation for the observed frequency peaks: K stands for the 
Keplerian frequency, L for the frequency of the lower branch of the Keplerian frequency, b 
for the break frequency. According to TO 99, the relation between the frequency of viscous 
oscillations Vy (symbol v in our table) and frequency i/b (Eq. 4) including equation (5) will 
serve as a tool for the identification of V), and v^. The last column of Table 2 contains Vh 
calculated according to formula (2) with VLq/2t{ = 380 Hz and also the angle S (calculated 
according to formula (5) for those cases when i/k and ul are measured simultaneously). 
The resulting values of S are shown in Figure 1. Indeed, the angle 6 shows little variation 
with uk and as expected is rather small 



5 = arcsin \{7r/no){iyl + ^^oA^)^^^('^i/^i^) 



(5) 



5 = 3.9° ± 0.2° 



(6) 
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With knowledge of 5 for the source 4U 1702-42, we attempted to interpret the remaining 
observed frequencies in Table 2. For July 26.60-26.93, uk and Ph are not present. Three 
frequencies 10.8, 32.5 and 80.1 Hz were observed on this day. If 32.5 Hz should be the ul 
then the second harmonic 2i/i = 65 Hz. But instead we have 80.1 Hz, which is significantly 
higher than the expected 21/^. On the other hand, from the classification Figure 4 presented 
in T099, we know that in the vicinity of i^k = 800 — 900 Hz, the two frequencies and z/y 
may come close to one another. Thus assuming that Uy — 32.5 Hz, we find that Uh — 10.8 Hz 
satisfies the theoretical relation (4) within the observational errors. We identify 2ul — 80.1 
Hz and also 21/^ = 85.5 Hz for July 21 and Uh = 12 Hz for the same day. From formula (3) 
we derive for both cases assuming the fixed S and fl/27: = 380 Hz. This identification 
holds best for 5 = 4.3° ± 0.1" which is the angle derived for — 769 Hz and = 40.1 Hz 
on July 30 (see the last fine of Table 2). The consistency of our identification is shown in 
Figure 2 where Uh, i'l, 2i'l and Ub are plotted as functions of (solid lines are theoretical). 
Figure 2 shows that as expected Au — Vh — decreases as vk increases. With the small 
number of observational points, the theoretical curves in Figure 2 are, in fact, a prediction 
of our model which can be viewed as a challenge for those observers who choose to extend 
the observational base for the source 4U 1702-42 to verify our model. It is possible that the 
difference between 5 = 3.8" and observed 5 = 4.3" is attributed to Q = const— assumption. 
Further measurement of and vk shall allow to reconstruct the Jl— profile and improve 
our predicative capability. 



4. Discussion and Conclusions 

Our model (OT99 and T099) reveals the physical nature of ~ 35 Hz oscillations in 
the source 4U 1702-42. As a low Keplerian branch of oscillations in the rotating frame 
of reference, this phenomenon should have an observational invariant - namely the angle 
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5. Within the first order approximation this angle can be viewed as a global parameter 
describing the inclination of the magnetospheric equator to the equatorial plane of the 
disk. Measured locally for different radial distances (therefore different uk), S may vary 
considerably unless indeed the observed oscillations correspond to the predicted low 
Keplerian branch. The constancy of 5 shown in Figure 1 allows us to interpret the ~ 35 
Hz oscillations as the low branch described by equation (3). Knowledge of the angle S is 
essential for the classification of the QPO resonances and the understanding of the physical 
nature of this phenomenon. The parameter 5 is critical in evaluating the differences between 
the spectra of different sources. For Sco X-1, oscillations with frequency ~ 45 Hz are 
identified as belonging to the ul branch (OT99). The higher observed frequencies of the 
for Sco X-1 are mainly a result of a larger angle 5 — 5.5° ± 0.5°. 

We believe that the angle 5 as a fundamental geometric parameter of the neutron star 
magnetosphere will be found eventually for other sources with kHz QPO. The source 4U 
1702-42 is the second source for which S has been inferred. 

The authors are grateful to J. Fainberg and R.G. Stone for discussions and suggestions. 
The comments of the referee which led to an essential improvement of the paper are 
appreciated. 
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VK (Hz) 


(Hz) 


f}/27r (Hz) 


657. 


1000. ± 7 


377. ± 5 


702. 


1038. ± 7 


382. ± 5 


769. 


1085. ± 11 


383. ± 8 



Table 1: Observed twin peak frequencies {pk and Uh) and calculated rotational frequency 
Q/27r for 4U 1702-42 (courtesy C. Markwardt). 
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Day of July, 1997 


Frequency (Hz) 


Classification 


Calculation 


19.37-19.86 


902. 


K 


Uh = 1179 Hz 




39. ± 3.4 


L 


S = 3.8° 


21.02-21.17 


85.6 ±3.7 


2L 






12. ± 1 


b 




26.60-26.93 


80.1 ±0.5 


2L 






32.5 ±1.2 


V 






10.8 ±0.4 


b 




30.46-30.72 


722. 


K 


Uh = 1048 Hz 




Au = 333. ± 5 




S = 3.8° 




34.7 ±0.6 


L 




30 


657. 


K 


Uh = 1005 Hz 




32.9 ±0.6 


L 


5 = 3.8° 


30 


702. 


K 


Uh = 1035 Hz 




34.8 ±0.6 


L 


S = 3.8° 


30 


769. 


K 


= 1081 Hz 




10.1 ±0.8 


L 


S = i.3" 



Table 2: Observed frequencies of 4U 1702-42 and possible theoretical interpretation. Angle 
5 is calculated from formula (5). Observations are cited from Markwardt et al. 1999 
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Fig. 1. — Inferred angle 5 between the rotational frequency fl and the normal to the plane 
of Keplerian oscillations as a function of observed Keplerian frequency for the source 4U 
1702-42. 



Fig. 2. — Classification of QPO in the atoll source 4U 1702-42. Sohd hues are theoretical 
curves. 
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